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FOREWORD 

The  project  discussed  in  this  technical  report  was  performed  under 
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results  of  analytical  investigation  into  the  behavior  of  plates  acoustically 
loaded  into  the  large  deflection  regime.  There  are  two  volumes  reporting 
this  research  effort:  Volume  I  consists  of  the  mathematical  formulations 
and  numerical  results;  Volume  II  contains  the  computer  codes.  The  study  was 
conducted  at  the  Department  of  Mechanical  Engineering  and  Mechanics,  Old 
Dominion  University,  Norfolk,  Virginia.  The  work  was  monitored  under  the 
supervision  of  Howard  F.  Wolfe,  AFWAL/FIBED,  Technical  Manager,  Acoustics 
and  Sonic  Fatigue  Group,  Structural  Integrity  Branch,  Structures  and 
Dynamics  Division,  Flight  Dynamics  Laboratory  of  the  Air  Force  Wright  Aero¬ 
nautical  Laboratories,  Wright-Patterson  AFB,  and  Dr.  Anthony  K.  Amos, 
AFOSR/NA,  Program  Administrator,  Directorate  of  Aerospace  Sciences,  Bolling 
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SECTION  I 
INTRODUCTION 

1.  BACKGROUND 

Acoustically  induced  fatigue  failures  in  aircraft  and  missile  operation 
have  been  a  design  consideration  for  over  25  years.  The  problem  was  intro¬ 
duced  with  the  advent  of  the  turbo  jet  engine  which  produced  high  intensity 
acoustic  pressure  fluctuations  on  aircraft  surfaces.  As  the  engine  perform¬ 
ance  requirements  increased,  the  intensity  of  the  acoustic  pressures  in¬ 
creased.  Airframe  minimun  weight  requirements  resulted  in  higher  stresses 
in  structural  components.  The  number  of  acoustic  fatigue  failures  began  to 
grow  at  a  rapid  rate  until  adequate  design  criteria  were  developed  and  used 
in  the  design  process. 

Similar  fatigue  failures  have  occurred  in  other  regions  of  high  inten¬ 
sity  pressure  fluctuations.  These  have  occurred  in  regions  of  separated 
flow,  behind  protuberances  into  the  flow  such  as  airbrakes,  and  close  to 
propeller  tips.  Failures  have  also  occurred  from  the  fluctuating  pressure 
induced  when  bomb  bay  doors  are  opened  during  high  speed  flight. 

The  oscillating  pressures  from  various  noise  sources  produced  a  reso¬ 
nant  response  of  the  structural  component  such  as  external  skin  panels, 
frames,  ribs  and  spars  which  results  in  rapid  stress  reversals  in  the  struc¬ 
ture.  If  these  stresses  have  sufficient  magnitude,  fatigue  failures  occur. 

Acoustic  fatigue  failures  have  resulted  in  unacceptable  maintenance  and 
inspection  burdens  associated  with  the  operation  of  the  aircraft.  In  some 
cases,  sonic  fatigue  failures  have  resulted  in  major  redesign  efforts  of 
aircraft  structural  components.  Therefore,  accurate  prediction  methods  are 
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needed  co  determine  the  acoustic  fatigue  life  of  structures.  Many  analyt¬ 
ical  and  experimental  programs  to  develop  sonic  fatigue  design  criteria, 
however,  have  repeatedly  shown  a  poor  comparison  between  analytical  and 
measured  RMS  maximum  strain  (References  1-9).  Deviations  in  excess  of  100 
percent  are  not  uncommon.  One  of  the  major  reasons  suspected  for  the 
discrepancy  was  that  the  panel  response  was  based  on  small  deflection  linear 
structural  theory,  whereas  the  test  panels  responded  with  large  deflections 
at  high  sound  pressure  levels  (SPL). 

Most  recently,  an  analytical  effort  (References  10,  11)  using  a  single¬ 
mode  approach,  and  a  test  program  (Reference  12)  have  demonstrated  that  the 
prediction  of  panel  random  response  is  greatly  improved  by  the  inclusion  of 
large  deflection  effects  in  the  analysis. 

Test  results  in  References  1,  2  and  12  also  have  shown  that  there  are 
more  than  one  mode  responding.  Multiple  modes  were  also  observed  by  Ublfe 
and  Wentz  (Reference  13)  and  White  (Reference  14)  in  experimental  investiga¬ 
tions  on  aluminum  fuselage  panels  and  carbon  fiber-reinforced  plastics 
(CFKP)  plates  under  acoustic  loading.  White  also  showed  that  the  fundamen¬ 
tal  mode  responded  significantly  and  contributed  more  than  one-half  of  the 
total  mean-square  strain  response;  higher  modes,  up  to  third  or  fourth 
modes,  account  for  95%  or  more  of  the  total  mean-square  strain  response.  In 
order  to  have  an  accurate  determination  of  the  random  response  of  a 
structure,  multiple  modes  should  be  used  in  the  analysis.  This  report.  Vol¬ 
ume  I,  presents  an  analytical  effort  on  large-amplitude  random  response  of  a 
clamped  rectangular  thin  plate  subjected  to  acoustic  excitation  using  multi¬ 
ple  (4)  modes  in  the  formulation.  The  developed  computer  software  is 
published  in  Volume  II  of  this  docunent. 
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2.  OBJECTIVES  OF  THE  INVESTIGATION 

The  primary  aim  of  this  research  is  to  investigate  the  large-amplitude, 
multimodal  response  of  clamped,  rectangular,  flat,  thin  panels  of  uniform 
thickness  subjected  to  broadband,  random  acoustic  excitation.  More 
specifically,  the  following  tasks  will  be  performed  in  accomplishing  this 
aim: 

(a)  Development  of  an  "quasi-exact"  solution  procedure  for  determining 
the  nonlinear  response  of  a  clamped  rectangular  plate  to  random  loading; 

(b)  Computerisation  of  the  solution  procedure  to  simplify  application 
to  the  nonlinear  random  response  problems  studied  in  this  invest igaf ion  and 
for  the  convenience  of  other  researchers  who  wish  to  use  and/or  extend  this 
work;  and 

(c)  Determination  of  the  RMS  maximum  panel  deflection,  the  KMS  maximum 
strain,  and  the  frequencies  of  vibration  of  an  acoustically  excitated  clamp¬ 
ed  plate  which  will  enable  the  practical  utilisation  for  eacimation  of  sonic 
fatigue  life. 
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SECTION  II 

MATHEMATICAL  FORMULATION  AND  SOLUTION  PROCEDURE 
I.  GOVERNING  EQUATIONS 

Tli*  governing  equations  of  a  thin  rectangular,  isotropic  plate  undergo¬ 
ing  large-deflection  notions,  neglecting  the  effects  of  both  inplane  and 
rotatory  inertia  forces,  are  (References  15,  16) 


L(w 


,F)  *  DV**w  •*■  phw,  +  gw, 

t  b  C 


h(F,  w,  +  F,  w,  -  2  F,  w,  ) 
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-  p(t)  -  0 
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where  a  comma  denotes  the  partial  differentiation  with  respect  to  the  corre¬ 
sponding  variable,  w  is  the  lateral  deflection,  D  is  the  flexural 
rigidity,  p  is  the  mass  density,  h  is  the  plate  thickness,  p  is  the 
pressure,  E  is  the  Young's  modulus,  and  g  is  the  viscous  damping.  The 
stress  function  F  is  defined  by 
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where  a0K>  °Qy  «nd  T0Xy  «r«  the  nembreae  stresses. 

Equations  (1)  and  (2)  together  with  a  complete  set  of  boundary  condi¬ 
tions,  define  the  problem.  These  equations  are  subjected  to  the  usual  out- 
of-plane  boundary  conditions  required  in  small  deflection  theory  (aero  edge 
deflection  and  slope  for  clamped  edges).  In  addition,  for  large  deflec¬ 
tions,  it  is  necessary  to  specify  inplane  boundary  conditions.  The  inplane 
boundary  conditions  must  be  specified  in  terms  of  a  combination  of  applied 
edge  load  and  displacements. 

The  particular  inplane  boundary  conditions  of  most  interest  in  the 
study  of  sonic  fatigue  of  structural  panels  is  the  one  in  which  the  edges 
are  restrained  from  movement,  that  is 


u  ■  0 

on  x  ■  0  and  a 

(6) 

v  *  0 

on  y  ■  0  and  b 

(7) 

where  the  plate  has  a  length  a,  width  b,  and  the  origin  of  the  coordin¬ 
ates  is  taken  at  one  corner  as  indicated  in  Figure  1.  The  numerical  results 
presented  in  the  next  section  are  all  based  on  these  inplane  boundary  condi¬ 
tions  . 

2.  SOLUTION 

Generally,  there  are  four  methods  which  sre  used  in  the  analysis  of 
mulci-degree-of-freehom  nonlinear  systems  under  random  excitations.  These 
are  the  Fokker-Planck  equation  approach,  the  perturbation  approach,  the 
finite  difference  approach  and  the  equivalent  (statistical  or  stochastic) 
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Figure  1.  Coordinates  and  Dimensions 


linearization  approach.  Application  of  these  methods  to  panels  of  complex 
configurations  were  discussed  in  Reference  4.  There  have  been  very  few 
attempts  at  solution  of  Equations  (1)  and  (2)  treating  plates  under  random 
loading. 

The  classical  approach  attempts  a  direct  solution  of  the  differential 
aquations,  thereby  allowing  a  reduction  in  numerical  approximations.  This 
approach  generally  involves  the  assumption  of  a  set  of  global  solution 
functions  which  are  simple  analytically,  usually  one  term  (References  10, 

11,  17  and  18)  and  satisfy  the  geometric  boundary  conditions.  The  differen¬ 
tial  equations  are  then  satisfied  approximately  by  use  of  a  Rayliegh-Ritz  or 
Galerkin  procedure.  Such  an  approach  allows  simple  solutions;  but  since  the 
shape  of  the  solution  functions  are  predetermined,  the  accuracy  of  the  pre¬ 
diction  is  limited  by  these  choices.  A  complete  description  of  the  deflec¬ 
tion  shape  will  not  be  able  to  emerge  since  the  shape  has  been  predeter¬ 
mined.  A  notable  exception  to  the  constrained  function  approach  was  the 
solution  procedure  developed  by  Levy  (Reference  19)  for  bending  of  a  simply 
supported  rectangular  plate  with  large  deflections. 

Levy's  procedure  involved  the  representation  of  the  stress  function  F 
and  the  lateral  deflection  function  w  by  double  infinite  Fourier  series. 
Such  a  representation  of  the  solution  functions  by  a  complete  set  of  "basic 
functions"  removes  the  major  limitation  of  constrained  solution  functions 
present  in  the  earlier  solutions.  The  relationship  between  the  coefficients 
of  the  stress  function  series  and  the  lateral  deflection  function  series  was 
determined  by  direct  substitution  of  the  two  series  into  the  compatibility 
equation  and  by  equating  coefficients  of  like  trigonometric  terms.  The 
equilibrium  equation  was  then  solved  by  the  same  procedure  of  substitution 


of  the  two  series  for  F  end  w  into  the  equilibrium  equation  and  equating 
coefficients  of  like  trigonometric  terms.  The  result  of  these  operations 
was  an  infinite  set  of  simultaneous  algebraic  equations  which  were  cubic  in 
the  deflection  function  coefficients.  This  set  of  equations  represents  a 
complete  solution  of  the  large  deflection  plate  bending  equations  for  the 
simply  supported  rectangular  plate.  Levy's  solution  is  referred  to  as  an 
"exact"  solution  by  Timoshenko  (Reference  15)  and  other  authors.  Of  course , 
the  exact  numerical  results  cannot  be  obtained  (this  would  require  solution 
of  an  infinite  number  of  simultaneous  equations) ,  but  since  the  series  is 
complete  and  the  solution  convergent,  the  degree  of  accuracy  of  the  results 
depends  only  on  the  patience  and  time  of  the  analyst.  Levy's  solution  pro¬ 
cedure  has  been  very  useful  to  the  research  community  as  evidenced  by  the 
many  investigators  which  have  made  use  of  and  extended  his  solution  to  study 
various  phenomena  associated  with  the  large  deflection  of  plates.  For 
example,  Paul  (Reference  20)  extended  Levy's  procedure  to  study  large  de¬ 
flection  behavior  of  clamped  plates  under  thermal  loading.  The  represent¬ 
ation  of  the  stress  function  and  the  plate  deflection  function  which  is  de¬ 
veloped  herein  will  also  follow  the  same  procedure  as  Levy  and  Paul. 

The  solution  will  be  written  in  terms  of  the  stress  function  F  and 
the  lateral  displacement  w,  with  the  relationship  between  these  variables 
determined  by  satisfying  the  compatibility  equation,  Equation  (2).  The 
lateral  deflection  w  is  represented  by  a  complete  set  of  functions  in  the 
form  of  a  "Fourier-type"  double  series. 

w(x,y,t)  ■  h  H  W  (t)  f  (x)  g  (y)  m,n  -  1,2,3...  (8) 

mn  m  n 

m  n 
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Tha  functions  fn(x)  and  gR(y)  are  chosen  to  be  linearly  independent  of  each 

other  and  such  that  the  geometric  boundary  conditions  of  sero  slope  and  de~ 

flection  along  the  edges  of  the  plate  are  satisfied  exactly.  The  W  (t) 

ran 

are  the  generalised  displacements  which  are  to  be  determined.  The  lateral 
boundary  conditions  are 


w  ■  w,  ■  0 
x 


on  x  ■  0  and  a 


w  ■  w,  ■  0 

y 


on  y  ■  0  and  b 


The  functions  f^Cx)  *nd  g^(y)  which  meet  the  above  requirements  are 


given  by: 


f  (x) 

m 


•nW 


(m+l)irx 


(n+l)in 

b 


To  help  perceive  the  numerical  behavior  for  the  component  functions  f  (x) 

fn 

and  8n^  which  make  up  the  series,  the  first  six  terms  of  the  series 

(Equation  8)  are  shown  in  Figures  2  through  7.  These  figures  are  plotted 

with  the  generalised  displacement  W  held  constant. 

mn 

Upon  examination  of  the  above  expression  of  the  lateral  deflection,  it 
is  found  that  the  compatibility  equation  (Equation  2)  can  be  identically 
satisfied  if  the  stress  function  F  is  taken  in  the  following  form: 
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Lateral  Deflection  Term  f1(x)g2(y) 


Figure  5.  Lateral  Deflection  Term  f2(x)g2(y) 
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P  x2 

—  |  •-  +  Eh2 


I  £ 

p  q 


Fpq  co8  008  p»q 


0,1,2...  (13) 


where  Che  constants  P  and  P  contribute  to  the  membrane  stress  and  will 

x  y 

be  discussed  in  more  detail  later.  Also,  it  should  be  noted  that  the  above 
form  of  the  stress  function  implies  zero  shear  stress  on  the  plate  bound¬ 
aries.  Direct  substitution  of  Equations  (8,  11,  12  and  13)  into  the 
compatibility  equation  Equation  (2),  performing  the  required  differentia¬ 
tions  and  multiplications,  and  equating  coefficients  of  like  trigonometric 

terms  yields  the  following  quadratic  relationship  between  F  and  W 

pq  mn 


F 


pq 


l 

(p2  /a  +  q2a)2 


HUB 
m  n  k  1 


pqrankl 


W  W.  . 
mn  kI 


(14) 


where  B  are  integers  and  a  ■  a/b.  Derivation  of  B  ,  ,  can  be 

pqmnkl  9  pqmnkl 

found  in  detail  in  Reference  20.  The  calculation  of  the  integers  Bpqnmjc^ 
is  accomplished  with  the  aid  of  a  computer  program,  a  listing  of  which  is 
presented  in  Volume  II  of  the  report.  A  tabulation  of  ®pqmn|lj  for  values 
of  the  subindices  p,  q,  m,  n,  k  and  1  from  zero  through  nine  has  been 
generated,  but  it  is  too  voluminous  to  be  included  in  the  report.  However, 
a  sample  of  this  tabulation  can  be  found  in  Reference  20.  Substitution  of 
Equation  (14)  into  Equation  (13)  yields  the  following  expression  for  the 
stress  function  in  terms  of  the  unknown  parameters  W  , 


P  y2  P  x2 

F  -  -  -j - ^  ♦  Eh2  £  £  [ 

P  q 


16 


*  •  e  v  UV  V  V  •*  »" .  •*%  <*  *’«*  »  *  <m  *  *V  \  s 


AAA 

A 


“"’■“SSEiift 


•  m*  k'*  *"■«•**  •  *  »  *  JS  *.»  * 


<■*. M<«  ra  *  J  WT  Tg^U.T  r  «t  V.^'CTVTOWVWIW"«r.-^-yi  ^.^TT^J^Wr.TT.  »%  IT. 


£li 


m 


—  ■  ■■  —  ■ — ——  I  E  E  I  B  , .  W  W. .  ]  cos  *--  cos  . 
(p2/a  +  q2«)2  m  n  k  1  P*"nkl  ran  kl  *  b 


Ptfl  “  0,1,2, .. . 

m,n,k,l  *  1,2,3, .. . 


(15) 


In  the  above  equation,  the  P  and  P  are  determined  by  applying  the  in- 

x  y 

plane  boundary  conditions,  and  the  Wffln  are  determined  by  solving  the  re¬ 
maining  governing  differential  Equation  (1). 

As  stated  earlier,  the  numerical  results  presented  in  this  study  are 
based  on  the  case  in  which  the  edges  are  restrained  from  movement  (Equa¬ 
tions  6  and  7).  In  order  to  apply  these  inplane  boundary  conditions, 
expressions  for  the  inplane  displacement  must  be  developed.  Combining  the 
membrane  strain-displacement  relations 


ox  'X 


u,  +  i-  v,2 


(16) 


*oy  "  v,y 


(17) 


Y  ■  u,  +  v,  +  w,  w, 
oxy  y  x  x  y 


(18) 


and  the  membrane  strain-stress  relations 
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and  re-arrangement  of  tha  raaulting  expressions  yiald 


/*  r^1  dx  ■  /*  [  4-  (F,  -  vF,  )  -  4  w,2]dx 

•  ox  \  1  B  *yy  *xx  2  *xJ 


(21) 


fb  §7  *y  "  lb  t  w  <F»w  -  VF»„J  -  7  w*i]dy 


.  -  E  '  *xx  *yy  2  "*yJ 


(22) 


where  E  is  the  Young's  modulus  and  v  is  the  Poisson's  ratio. 

Performing  the  indicated  differentiation,  integration  and  applying  the 
in-plane  boundary  condition  (Equations  6,  7)  to  each  of  the  above  equations 
respectively  yield: 


Pa  vp  a 


E  E 


X _ - 


Ix(W«n) 


(23) 


0  -  - 


P  b  vp  b 
-2~-  +  — 2L_ 
E  E 


I  (W  ) 
y  mn 


(24) 


The  terms  I  (W  )  and  I  (W  )  in  the  above  equations  are  defined  by 
x  mn  y  ran 


2  2 

1  (v  )  -  e  z  w  z  (w  ) 

x  ran  8a  mn  x  mn 

m  n 


(25) 


h2ir2 

1  (W  )  -  E  E  W  Z  (W  ) 
y  ran  8b  ran  y  mn 

ra  n 


(26) 
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(27) 


*_<*>  *  t(n+1)2  +  <n-l)2]W.  „  *  (n+1)2  W, 

y  ®n  2-o,n  2-o,n 

-(n-l)2W_  ,  +  2[(n+l)2  +  (n-l)2]W 

Z-o,n-Z  an 

-2(n+l)2W  -  2(n-l)2W  , 

o,n-z  a,n-2 

-[ (n+1)2  ♦  (n-i)2 ]W  +  (n+l)2W  , 

a-z,n  ra-z,n+z 

+(n-l)2W  -  n  ,  -  [ (n+1)2  +  (n-1)2 JW  « 

o-2, n-2  m+2, n 

+(n+l)2W  +  (n-l)2W  _ 

o+2, n+2  o+2, n-2 


(28) 
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where  W  •  0  for  ra  or  n  <  1. 

mn 


Solving  Equation*  (23)  and  (24)  ilmultaneoualy  yields  the  following 

expressions  for  the  average  edge  loads  P  and  P  for  the  inplane 

x  y 

boundary  conditions  in  which  the  edges  are  restrained  from  movement 


m 


'x  •  •  -pf*-  (b  *y  *  r  V 


rir (V.  *  r  Iy> 


(29) 


(30) 


With  the  lateral  deflection  given  by  Equation  (8)  and  the  stress 
function  given  by  Equation  (15),  the  equilibrium  equation,  Equation  (1),  is 
then  satisfied  by  applying  the  Bubnov-Galerkin  method  (Reference  21) 


/b  /*  L  (w,F)  ffga  dxdy  -  0  r,s  *  1,2,3... 


(31) 


Upon  substitution  of  L(w,F)  in  Equation  (31),  the  equations  are  written 


as 


C  C  +  »*>t  *  *  2"-*xyy  *  w-yyyy) 


M 


-h(F,  w,  +  F,  w,  -  2F,  w,  )-p]f  g  dxdy  ■  0 
yy  xx  ’xx  yy  ’xy  xy  r  a 


where  r,s  ■  1,2,3,... 


(32) 


Utilising  the  expressions  for  w,  f  ,  g  and  F  (Equations  8,  11,  12 

4  8 
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[K]  (W)  ♦  tC]  {W}  +  [K]l  {W}  ♦  {<5(W»  -  { p( c )}  (33) 

where  [Mj,  [C]  end  (K]^  ere  Che  generalized  mass,  damping  and  linear 
stiffness  matrices,  respectively,  and  ($(W)}  is  a  vector  function  which  is 
cubic  in  the  generalised  displacements  (W).  The  integrations  of  each  of 
the  terms  in  Equation  (32)  are  given  explicitly  in  Appendix  A.  The  mass  and 
linear  stiffness  matrices  can  be  easily  generated  from  equations  (Al)  to 
(A4) . 

If  the  acoustic  pressure  loading  p(t)  is  stationary  Gaussian,  is 
ergodic  and  has  a  zero  mean.  Application  of  the  Krylov-Bogoliubov-Caughey 
equivalent  linearization  method  (References  22-28)  yields  an  equivalent  set 
of  linearised  equations  to  Equation  (33)  as 

[M]  {W}  +  [C]  {W }  +  ((K]l  ♦  (K]gL)  {W>  -  (p(t)}  (34) 

or 

[Ml  (W)  +  [C]  (W)  +  IK]  {W}  -  (p(t)}  (35) 

where  [KK  is  the  generalized  equivalent  linear  stiffness  matrix  and 
[K]  ■  [Kj^  +  The  elements  of  are  evaluated  from  the  equation 

(Reference  23) 
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f^Jt.ij  ■  {  tsw^j  1  (36) 

and  thay  art  darivad  and  given  axplicitly  in  Appendix  B,  and  CM  stands 
for  tha  oparator  of  mathematical  axpactation.  Tha  approximate  generalised 
displacements  {W},  computed  from  tha  linearised  equation,  Equation  (35), 
are  also  Gaussian  and  approach  stationary  because  the  panel  motion  is 
stable. 

To  determine  the  mean-square  generalised  displacements  W2^  in  Equa¬ 
tion  (34),  an  iterative  process  is  introduced.  The  undamped  linear  equation 
of  Equation  (34)  is  solved  first,  which  requires  the  determination  of  the 
eigenvalues  and  eigenvectors  of  the  undamped  linear  equation 


•?  [M]  4>.  -  (KL4). 
J  J  L  J 


(37) 


where  ok  is  the  linear  frequency  of  vibration,  and  4K  is  the  normal 
mode  shape. 

Apply  a  coordinate  transformation,  from  the  generalised  displacements 
to  the  normal  coordinates  (4  modes  will  be  used  in  the  analysis,  n  "  4),  by 


{ W)  ■  4]  { q}  n<m 

mxl  raxn  nxl  (38) 


where  each  column  of  4]  is  a  normal  mode  4K. 
equation  of  Equation  (34)  becomes 


The  damped  linear 
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vhera  f* MoJ  ■  for  [M]  {♦] 


hc>)  -  WP  [c]  [*)  •  i  r*c».j  hM^ 


Is  k *JL  -  [♦J1  (K]L  (♦)  -  (Na2^  h*K) 


{p}  •  l*r  {p} 


Equation  (39)  is  uncoupled  in  the  normal  coordinates,  the  jth  row  of  Equa¬ 
tion  (39)  is 


’j  * 2  *  i"  ■  sf 


The  mean-square  response  of  modal  amplitude  is 


q2.  -  Jm  Sp(w)  |H.(u)  j2  dw 

J  ®  "  J 


where  Sp(m)  is  the  spectral  density  function  of  the  excitation  Pj(t),  and 
the  frequency  response  function  is  given  by 


"  M.(fttf  -  w2  ♦  2iC  .w .“) 
J  J  J  J  J 
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For  lightly  damped  (C  »  c/c  <  0.03)  structural  ths  response  curves  will 
be  highly  peeked  et  Wj .  the  integration  of  Equation  (43)  can  be  greatly 
simplified  if  the  forcing  spectral  density  function  Sp(w)  can  be  consid¬ 
ered  to  be  constant  in  the  frequency  band  surrounding  the  linear  resonance 
peak  sij  |  so  that 


— ■  tSttOo.) 


<6  « 


3 


(47) 


The  covariance  matrix  of  the  linear  generalised  displacements  is 


S„(<*».) 


(48) 


The  diagonal  terms  of  Iwmn\i^L  *re  the  raean“8<luare  linear  generalized 

displacements  W®n<  This  initial  estimate  of  expected  value  on  generalized 
displacements  now  can  be  used  to  compute  the  generalized  equivalent  linear 
stiffness  matrix  £ K] through  Equation  (36)  or  Equations  (B4  to  B7).  The 
undamped  linearized  equation  of  Equation  (34)  is  solved  again 


1)2  (HlWj  -  <[K)l  ♦  [«,,_)  (♦). 


(49) 


where  is  the  equivalent  linear  or  nonlinear  frequency  of  vibration,  and 
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{♦> j  is  tha  aaaoeiaead  equivalent  linear  normal  mode.  Then  Equation  (35) 
ie  transformed  again  to  the  normal  coordinates  and  it  has  the  form  as 

MN  (q)  ♦  r-Cg  (1)  ♦  MU]  (q)  -  (P(t))  (50) 

in  which 

MN  -  i*]1  <[K]L  4  IK]EL)W]  -  MN  (5i) 

The  jth  row  of  Equation  (30)  is 

.  P. 

'ij  ♦  2Cj  Vj  *  "fo  ■  (52> 

and  the  displacement  covariance  matrix  is 


[W  W.  ,  J 
ran  klJ 


I  (♦)  .  J — •« 

3  J 


Sp(fl.)  T 

— — J — (♦>* 
m2  r  M  o2  J 


The  diagonal  terms  of  (W  W. .  ]  are  the  mean-square  generalized  displace- 

ran  ki 


ments  W2 ^ .  As  the  iterative  process  converges  on  the  iter-th  cycle,  the 
mn 

relations 


(q*)..  B  (q*)...  . 

j  iter  j  iter-1 
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become  satisfied.  In  the  numerical  raaulcc  presented  in  tha  following 
■action!  convarganca  it  considered  achiavad  whan  tha  diffaranca  of  tha  RMS 
diaplacamanta  aatiafiea  tha  ralation 


<BMS  Vi,.,  *  <"18  Vi«.,-i  t 

- 10  f0r  U 


Once  the  RMS  ganaralisad  diaplacamanta  are  datermined,  the  RMS  deflec¬ 
tion  of  the  panel  and  the  RMS  maximum  strain  can  be  determined  from  Equa¬ 
tions  (8,  19,  20)  and  the  bending  strain-displacement  relations  as  presented 
in  the  next  section. 

In  practice,  the  spectral  density  is  usually  given  or  measured  in  terms 
of  the  frequency  f  in  Harts.  To  convert  the  spectral  density  from  the 
analytical  Sp(w)  to  practical  Sp(f),  the  following  relations  are  used 
(Reference  29) 


u  -  2"f 


8p(w)  -  Sp(f)/2ff 


The  pressure  spectral  density  function  Sp(f)/2ir  then  has  the  units 
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(Pa)2 /Ha  or  (p«i)2/Ha  and  tha  diaplaeamant  covariance  matrix,  Equation 
(53),  becoraea 


W*t‘  <♦> •  r~~-  <♦>: 


“  'T*  1  - - —  iV/  1 

3  j  "j'i'oj'j  J 


(60) 


whore  S ^  ia  a  nondlmenaional  forcing  apectral  denaity  parameter  defined  as 


and 


8p(f) 


p2h4(0/pb4)s/2 


w2  .  \2  t 

o  ^phb1^ 


(61) 


(62) 


1)2  "  *2  (s5Pr) 


(63) 


M/ph2 


(64) 


The  frequency  parameters  \q  and  A,  and  the  mass  parameter  M  are  all 
dimensionless  quantities.  This  nondioansional  spectral  density  parameter  of 
acoustic  pressure  excitation  will  be  used  for  the  numerical  results 
presented  in  the  next  section. 


3. 


DEFLECTION  AND  STRAINS 

From  Equation  (8),  the  mean-square  deflection  of  the  plate  is  given  by 


mmmmm 
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(65) 


Where  the  averages  W  W,  -  ■  5{W  W.  . ]  have  been  determined  in  Equation 

mn  Ki  .  mn  ki 

(53).  The  mean-square  center  deflection  can  be  obtained  from  Equation  (65). 

For  a  clamped  rectangular  plate  (a  >  b),  the  maximum  bending  strain 
occurs  at  the  extreme-fiber  (z  ■  ±  h/2)  at  the  midpoint  (x  ■  a/2)  of  the 
long  edge  and  perpendicular  (e^)  to  the  long  edge.  The  bending  strain  at 
the  extreme-fiber  and  in  the  y-direction  is 


e  -  ±  H  §!* 

*  23y* 


(66) 


or 
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(67) 


The  membrane  strain  in  the  y-direction,  Equations  (15,  20,  29  and  30), 
in  terms  of  generalized  displacements  Wmn  is 
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m  n  k  1  p  q 


(68) 


where 


pq 


1>2(Vq2  -  p2/q2)  ^  £[x  ^ 

e  b 


(p2/o  +  q2  a)2 


(69) 


The  total  strain  is  then  given  by 


e  -  e  +  e. 
y  oy  by 


(70) 


and  the  mean-square  strain  becomes 


y  oy 


e2  +  2C  [«  e.  ]  +  e 

Aif  •»  a  \r  hwJ 


oy  byJ  by 


(71) 


For  Gaussian  random  processes  with  zero-mean,  the  generalized  displacement 
follow  the  relations 


cfw.  .W.  .W  ]  -  0 

1  1J  kl  mnJ 


(72) 


cfw.  ,W_W  W  1  -Cfw,  .W..1  c  [w  W  ] 
1  ij  kl  mn  rsJ  1  ij  klJ  1  inn  rsJ 


♦  cfw.  .W  ]  c  [w..W  1  +  c[w..w  ]  cfw-.w  ]  (73) 

1  ij  mnJ  *■  kl  rsJ  1  ij  rsJ  1  kl  mnJ 


and  the  RMS  maximum  strain  can  be  determined  from  Equation  (71). 
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SECTION  III 

RESULTS  AND  DISCUSSION 

1.  GENERAL  COMMENTS 

Since  the  solution  developed  in  Section  II  is  of  the  infinite  series 
type,  it  is  important  to  understand  its  convergence  behavior.  The  experi¬ 
ence  gained  in  Reference  20  on  convergence  of  the  solution  for  the  large 
deflection  under  thermal  loading  aids  in  truncation  of  the  series  in  a 
rational  manner. 

The  convergence  of  the  solution  is  examined  through  a  study  of  square 
and  rectangular  panels  subjected  to  a  normal  incidence  acoustic  impingement. 
Since  the  loading  is  symmetric,  therefore,  only  symmetric  generalized  dis¬ 
placements  are  retained  in  the  transverse  deflection  function.  The  partic¬ 
ular  generalized  displacements  that  were  chosen  to  be  nonzero  in  the  con¬ 
vergence  study  are  shown  in  Table  1.  Numerical  results  are  provided  for  the 
mean-square  plate  center  deflection  and  the  mean-square  maximum  strain  for  a 
lateral  deflection  series  that  varies  in  size  from  four,  six,  ten  and  to 
fifteen  terms.  The  maximum  strain  occurs  at  the  extreme-fiber  (z  ■  ±  h/2) 
perpendicular  to  the  long  edge  (e^)  and  at  the  midpoint  (x  *  a/2).  In 
the  results  presented,  the  white  noise  excitation  is  band-limited  with  a 
frequency  bandwidth  of  25  Hz  to  6,000  Hz  as  shown  in  Figure  8,  the  damping 
ratio  is  assumed  to  be  constant  for  all  four  normal  modes  and  Poisson's 
ratio  is  equal  to  0.33.  The  computer  program  is,  however,  much  more  flex¬ 
ible  that  the  spectral  density  of  pressure  excitation  could  be  varying  slow¬ 
ly  in  the  neighborhood  of  equivalent  linear  frequency  0j  and  also  the 
damping  ratios  could  be  different  for  different  normal  coordinates. 
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TABLE  1.  GENERALIZED  DISPLACEMENTS  FOR  CONVERGENCE  STUDIES. 


Number 

0 

Hi 

IT 

• 

H 

5 

* 

Generalised 

Displacements 

4 

6 

10 

_15 

W11 

X 

X 

X 

X 

w13 

X 

X 

X 

X 

»31 

X 

X 

X 

X 

«33 

X 

X 

X 

X 

W1S 

X 

X 

X 

^1 

X 

X 

X 

W35 

X 

X 

*%3 

X 

X 

*17 

X 

X 

*71 

X 

X 

*55 

X 

*17 

X 

*73 

X 

*19 

X 

w91 

X 

2.  CONVERGENCE  OP  THE  NONLINEAR  RANDOM  RESPONSE 

The  convergence  of  the  solution  technique  was  examined  to  determine  the 
degree  of  accuracy  possible  with  a  highly  trancated  transverse  deflection 
function  series,  the  mean-square  nondimensional  center  deflection  versus 
the  nondimensional  spectral  density  paraneter  Sf  using  4,  6,  10  and  15 
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Figure  8.  Band— Limited  White  Noise  Excitation 


terra#  in  the  deflection  function  foe  a  square  end  e  rectangular  (a  ■  2) 
plate  are  shown  in  Figure*  9  and  10,  respectively.  It  clearly  indicate* 
that  a  6-terra  solution  gives  accurate  results  for  the  nonlinear  maximum  de¬ 
flection  while  a  4-term  solution  will  provide  very  accurate  linear  results. 
Figures  11  and  12  show  the  maximum  mean-square  strain  versus  the  nondimen- 
sional  pressure  spectral  density  for  the  square  and  rectangular  (a  ■  2) 
plates,  respectively,  using  4,  6,  10  and  15  terms  in  the  deflection 
function.  It  can  be  seen  that  the  convergence  of  the  mean-square  strain  is 
much  slower  as  compared  with  that  of  the  mean-square  deflection.  More 
number  of  terms,  therefore,  are  needed  in  the  deflection  function  for  accur¬ 
ate  determination  of  strains.  Results  of  maximum  mean-square  strain  based 
on  small  deflection  linear  theory  are  also  given  in  the  figures.  The  use  of 
linear  theory  would  lead  to  poor  estimation  of  panel  fatigue  life. 

3.  EFFECTS  OF  PANEL  DAMPING 

Figures  13  and  14  show  the  mean-square  nondimens ional  center  deflection 
versus  the  nondimens ional  spectral  density  of  acoustic  pressure  excitation 
for  rectangular  panels  of  length-to-width  ratio  o  ■  l  and  2  with  the  damp¬ 
ing  ratio  equal  to  0.009,  0.018  and  0.027.  Figures  15  and  16  show  the 
maximum  mean-square  nondimens ional  strain  versus  the  nondiraens ional  pressure 
spectral  density  for  rectangular  plates  of  aspect  ratios  of  1  and  2  with 
different  panel  damping  ratios.  Ten  terras  were  included  in  the  deflection 
function  to  generate  those  results  shown  in  Figures  13  through  16.  It  is 
clear  from  the  figures  that  the  precise  determination  of  panel  damping  is 
very  important. 
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SECTION  IV 
CONCLUDING  REMARKS 

A  "quasi-exact"  analytical  procedure  was  developed  £or  the  prediction 
of  the  nonlinear  random  response  of  rigidly  clamped  rectangular  panels  sub¬ 
jected  to  broadband  noise  excitation.  Multiple  (4)  modes  were  employed  in 
the  formulation.  A  computer  program  has  been  developed  to  aid  in  the 
determination  of  RMS  center  deflection,  RMS  maximum  strain  and  equivalent 
linear  frequency  at  given  pressure  spectral  density  of  excitation. 

The  convergence  of  the  present  analytical  method  was  demonstrated 
through  a  detailed  study  of  a  square  and  a  rectangular  (a  ■  2)  plate  under 
normal  incidence  acoustic  impingements.  Results  were  presented  for  both 
panels  with  the  truncated  lateral  deflection  function  series  of  4,  6,  10  and 
IS  terms .  This  study  revealed  that  convergence  is  very  rapid  for  the  deter¬ 
mination  of  central  deflection,  but  much  slower  in  the  prediction  of 
strains.  Accurate  mean-square  deflections  can  be  obtained  with  the  use  of 
six  terms  in  the  deflection  function,  while  it  is  necessary  to  consider  as 
many  as  10  or  more  terms  for  the  accurate  determination  of  the  strains. 

This  computed  RMS  strain  and  equivalent  linear  frequency,  in  conjunction 
with  strain  versus  cycles  to  failure  (S-N)  data,  should  be  used  for  the 
etimation  of  service  life. 

In  the  numerical  examples  presented,  a  constant  damping  ratio  was  used 
for  all  excitation  levels.  The  total  damping  includes  the  acoustic  radi¬ 
ation  damping,  panel  edge  damping  and  the  material  damping  in  the  panel 
itself.  Nonlinear  damping  phenomena,  however,  have  been  observed  in  experi 
ments.  Strain  response  data  for  an  aluninun  panel  of  3.94  cm  x  3.94  cm  x 
0.23  mm  at  three  different  overall  sound  pressure  levels  (130,  142  and  160 
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dB)  obtained  in  Reference  12  are  shown  in  Figure  17.  At  the  low-excitation 
levels,  the  modal  responses  of  the  panel  can  be  clearly  seen.  At  the  high¬ 
est  level,  broadening  of  the  response  peak  due  to  nonlinear  damping  is 
observed.  Therefore,  more  effort  is  needed  to  better  understand  the  effects 
of  nonlinear  damping  on  panel  response. 
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APPENDIX  E 


GENERALIZED  EQUIVALENT  LINEAR  STIFFNESS  MATRIX  [K]bl 
The  • laments  of  ths  equivalent  linssr  stiffnsss  matrix  are  derived  from 
the  expression 


(Bl) 


The  function  is  associated  with  the  integrals  after  applying  the 

Bubnov-Galerkin  approach  to  the  aquation  of  motion  in  deflection,  and  it  can 
be  expressed  as  the  sum  of  five  subfunctions  as 


*r.  '  ‘S'  C  1‘  <r>y,"*«x  *  f’xx“'yy  ‘  2P•xy“•xy)fr*.dx'll, 


‘  »„>1  *  <»„)2  ♦  (Sr.>3  ♦  <«„>5 
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The  equivalent  linear  stiffness  matrix  will  consist  of  five  submatrices 

as 
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and  [wmn  Wkl]  i«  the  covariance  matrix  of  the  generalised  displacements. 

rha  terms  of  [w  W. .]  ere  the  maan~square  displacements  W2  . 

“*u  A  mn 

Expression  for  which  is  similar  to  has  also  been 

derived. 


REFERENCES 


1.  J.  Soovera,  "The  Effect  of  Acoustic-Thermal  Environments  on  Advanced 
Composite  Fuselage  Panels 24th  Structures,  Structural  Dynamics  and 
Materials  Conference,  Lake  Tahoe,  NV,  May  1983,  pp.  466-472. 

2.  J.  Soovera,  "8onic  Fatigue  Tasting  of  an  Advanced  Composite  Aileron," 
Journal  of  Aircraft,  Vol.  19,  April  1982. 

3.  I.  Ho l chouse,  "Sonic  Fatigue  Design  Techniques  For  Advanced  Composite 
Aircraft  Structures,"  AFWAL-TR-80-3019,  AF  Wright  Aeronautical  Labora¬ 
tories,  WPAFB,  OH,  April  1980. 

4.  C.  Mei,  "Large  Amplitude  Response  of  Complex  Structures  Due  to  High 
Intensity  Noise,"  AFFDL-TR-79-3028,  AF  Flight  Dynamics  Laboratory, 
WPAFB,  OH,  April  1979. 

5.  M.  J.  Jacobson,  "Sonic  Fatigue  Design  Data  for  Bonded  Aluminum  Aircraft 
Structures,"  AFFDL-TR-77-45,  WPFAB,  OH,  1977. 

6.  R.  C.  W.  Van  der  Heyde,  and  D.  L.  Smith,  "Sonic  Fatigue  Resistance  of 
Skin-Stringer  Panels,"  AFFDL-TM-73-149-FYA,  WPAFB,  OH,  1974. 

7.  M.  J.  Jacobson,  "Advanced  Composite  Joints;  Design  and  Acoustic  Fatigue 
Characteristics,"  AFFDL-TR-71-126,  WPAFB,  OH,  1972. 

8.  L.  D.  Jacobs,  and  D.  R.  Lagerquist,  "Finite  Element  Analysis  of  Complex 
Panel  to  Random  Loads,"  AFFDL-TR-68-44,  WPAFB,  OH,  1968. 

9.  G.  E.  Fitch,  T.  R.  Dutko,  L.  M.  Brennan,  A.  G.  Tipton,  P.  M.  Belcher, 

P.  Wang,  and  P.  A.  Clawson,  "Establishment  of  the  Approach  to,  and 
Development  of  Interim  Design  Criteria  for  Sonic  Fatigue,"  ASD-TDR- 
62-26,  WPAFB,  OH,  1962. 

10.  K.  R.  Wentz,  D.  B.  Paul,  and  C.  Mei,  "Large  Amplitude  Random  Response 
of  Symmetric  Laminated  Composite  Plates,  "Shock  and  Vibration  Bulletin 
52,  May  1982,  pp.  99-111. 

11.  C.  Mei,  "Response  of  Nonlinear  Structural  Panels  Subjected  to  High  In¬ 
tensity  Noise,"  AFWAL-TR-8Q-3018,  AF  Wright  Aeronautical  Laboratories, 
WPAFB,  OH,  March  1980. 

12.  C.  Mei,  and  K.  R.  Wentz,  "Analytical  and  Experimental  Nonlinear  Response 
of  Rectangular  Panels  to  Acoustic  Excitation,"  23rd  Structures, 
Structural  Dynamics  and  Materials  Conference,  New  Orleans,  LA,  May  1982, 
pp.  514-520. 

13.  H.  F.  Wolfe,  and  K.  R.  Wentz,  "Acoustic  Fatigue  Test  Evaluation  of  Ad¬ 
hesively  Bonded  Aluminun  Fuselage  Panels  using  FM73/BR126  Adhesive/ 
Primer  System,"  24th  Structures,  Structural  Dynamics  and  Materials 
Conference,  Lake  Tahoe,  NV,  May  1983,  pp.  626-631. 


14.  R.  G.  White,  "Comperieon  of  the  Statistical  Properties  of  the  Aluminum 
Alloy  and  CPRP  Plates  to  Acoustic  Excitation,"  J.  Composites,  Oct.  1978, 
pp.  251-258. 

15.  S.  Timoshenko,  and  S.  Woinowsky-Krieger,  "Theory  of  Plates  and  Shells," 
McGraw-Hill,  415-428,  1959. 

16.  H.  N.  Chu,  and  G.  Herrmann,  "Influence  of  Large  Amplitudes  on  Free 
Flexural  Vibrations  of  Rectangular  Elastic  Plates,"  J.  Applied  Mechan¬ 
ics,  Vol.  23,  December  1956,  pp.  532-540. 

17.  R.  E.  Herbert,  "Random  Vibrations  of  Plates  with  Large  Amplitude,"  J. 

Applied  Mechanics,  Vol.  32,  Sept.  1965,  pp.  547-552.  , 

18.  Y.  K.  Lin,  "Response  of  a  Nonlinear  Flat  Panel  to  Periodic  and  Random¬ 
ly-Varying  Loadings,"  J.  Aerospace  Sci.,  Sept.  1962,  pp.  1029-1033, 

p.  1066. 

19.  S.  Levy,  "Bending  of  Rectangular  Plates  with  Large  Deflections,"  NACA 
Report  737,  1942. 

20.  D.  B.  Paul,  "Large  Deflections  of  Clamped  Rectangular  Plates  with 
Arbitrary  Temperature  Distributions,"  AFWAL-TR-81-3003,  Volume  I,  AF 
Wright  Aeronautical  Laboratories,  WPAFB,  OH,  Feb.  1982. 

21.  A.  S.  Volmir,  "Flexible  Plates  and  Shells,"  AFFDL-TR-66-127,  AF  Flight 
Dynamics  Laboratory,  WPAFB,  OH,  April  1967. 

22.  P.  T.  D.  Spanos,  and  W.  D.  Iwan,  "On  the  Existence  and  Uniqueness  of 
Solutions  Generated  by  Equivalent  Linearization,"  Int.  J.  Nonlinear 
Mech.,  Vol.  13,  1978,  pp.  71-78. 

23.  T.  S.  Atalik,  and  S.  Utku,  "Stochastic  Linearization  of  Multi-Degree-of 
Freedom  Nonlinear  Systems,"  Earthquake  Eng.  Struct.  Dyn.,  Vol.  4,  1976, 
pp.  411-420. 

24.  W.  D.  Iwan,  and  I.  M.  Yang,  "Application  of  Statistical  Linearization 
Technique  to  Nonlinear  Multidegree-of-Freedom  Systems,"  J.  Applied 
Mechanics,  Vol.  39,  June  1972,  pp.  545-550. 

25.  T.  K.  Caughey,  "Nonlinear  Theory  of  Random  Vibrations,"  in  Advances  in 
Applied  Mechanics,  Edited  by  C.  S.  Yih,  Vol.  11,  Academic  Press,  1971, 
pp.  209-253. 

26.  E.  T.  Foster,  Jr.,  "Semilinear  Random  Vibrations  in  Discrete  Systems," 

J.  Applied  Mechanics,  Vol.  35,  Sept.  1968,  pp.  560-564. 

27.  T.  K.  Caughey,  "Equivalent  Linearisation  Techniques,"  JASA,  Vol.  35, 

1963,  pp.  1706-1711. 

28.  N.  Krylov,  and  N.  Bogoliubov,  "Introduction  to  Nonlinear  Mechanics," 
translated  by  Lefshletz,  S.  in  Annals  of  Mathematical  Studies,  No.  11, 
Princeton  University  Press,  N.  J. ,  1943. 


56 


29.  F.  F.  Rudder,  Jr.  end  H.  £.  Plumblee,  Jr.,  "Sonic  Fatigue  Design  Guide 
for  Military  Aircraft,"  AFFDL-TR-74-112,  AF  Flight  Dynamics  Laboratory, 
WPAFB,  OH,  May  1975,  p.  265-266. 


a 


"if  U  S  GOVERNMENT  PRINTING  OFFICE-7S9-O02/727 

57 


